*********************************************
* Forced migration and food crises		
* 										
*		Carril-Caccia, Paniagua and Suarez-Varela	
*										
* 		JOPE 2025								
*********************************************

*************************************************************************************************************
*************************************************************************************************************
*************************************************************************************************************
********************* Organization **************************************************************************
*************************************************************************************************************
*************************************************************************************************************
*************************************************************************************************************

**** The following code replicates the code for of the estimates presented in the paper in the following order:

* Table 1
* Figure 5
* Table 2
* Figure 6
* Table 3 (Tables A8, A9 & A10)
* Table A4
* Table A5
* Table A6
* Table A7
* Table 4

clear
clear matrix
set more off

cd "C:\Data"

use replication_migfood.dta, clear

global controlsfull intfoodcrisis intany_crisis intgdppcgr  intL1VoiceandAccountability intL1death_pop  rta lnmigration  INT lngdppc_s lngdppc_d  lnpop_s lnpop_d lndist contig comlang_off  colony 
global control1 intfoodcrisis
global control2 intfoodcrisis rta lnmigration 
global control3 intfoodcrisis intany_crisis rta lnmigration intgdppcgr
global control4 intfoodcrisis rta lnmigration intany_crisis intgdppcgr intL1VoiceandAccountability intL1death_pop
global controltab intfoodcrisis intany_crisis intgdppcgr  intL1VoiceandAccountability intL1death_pop rta lnmigration


global controlnf intany_crisis intgdppcgr  intL1VoiceandAccountability intL1death_pop rta lnmigration
global control4nc intany_crisis intgdppcgr  intL1VoiceandAccountability intL1death_pop rta lnmigration

global control1q quarterscount
global control2q quarterscount rta lnmigration 
global control3q quarterscount intany_crisis rta lnmigration intgdppcgr
global control4q quarterscount rta lnmigration intany_crisis intgdppcgr intL1VoiceandAccountability intL1death_pop
global controltabq quarterscount intany_crisis intgdppcgr  intL1VoiceandAccountability intL1death_pop rta lnmigration


global control4ncnoint  any_crisis gdppcgr  L1VoiceandAccountability L1death_pop 
global control4nc2 intany_crisis intgdppcgr  intL1VoiceandAccountability intL1death_pop
global bilat INT lndist contig comlang_off colony 
global country  lngdppc_s lngdppc_d lnpop_s lnpop_d 
global bilatt rta lnmigration


global controltab2  intany_crisis_t2  intgdppcgr_t2  intL2VoiceandAccountability   intL2death_pop  L2rta lnmigration
global controlnf intany_crisis intgdppcgr  intL1VoiceandAccountability intL1death_pop rta lnmigration
global controllag1 L2rta lnmigration
global controlag2 intany_crisis_t2  intgdppcgr_t2  L2rta lnmigration
global controlnf intany_crisis intgdppcgr  intL1VoiceandAccountability intL1death_pop rta lnmigration





cd "C:\Results"


*************************************************************************************************************************************************************
*************************************************************************************************************************************************************
*********************************************************Table 1*********************************************************************************************
*************************************************************************************************************************************************************
*************************************************************************************************************************************************************
local mylistd     application_IDP
foreach d of local mylistd {

eststo, ti("`d'"): ppmlhdfe application_IDP  $control1  INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store `d'1



eststo, ti("`d'"): ppmlhdfe application_IDP  $control2 INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store `d'2



eststo, ti("`d'"): ppmlhdfe application_IDP  $control3  INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store `d'3


eststo, ti("`d'"): ppmlhdfe application_IDP  $control4  INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store `d'4

}



esttab application_IDP1 application_IDP2 application_IDP3 application_IDP4	 using "Table1.tex", pr2 ///
replace nocons star(* 0.10 ** 0.05 *** 0.01) b(3) se(3)  ///
keep($controltab)  ///
order($controltab)   ///
label nonotes s(num_n FE3 FE5 FE6 FE7 IDP fo xt, label("Observations" "Country pair FE" "Origin-year FE" "Destination-year FE" "Border-year FE" "IDP" "foodt-1 ties INT" "xit-1 times int"))

esttab application_IDP1 application_IDP2 application_IDP3 application_IDP4	 using "Table1.rtf", pr2 ///
replace nocons star(* 0.10 ** 0.05 *** 0.01) b(3) se(3)  ///
keep($controltab)  ///
order($controltab)   ///
label nonotes s(num_n FE3 FE5 FE6 FE7 IDP fo xt, label("Observations" "Country pair FE" "Origin-year FE" "Destination-year FE" "Border-year FE" "IDP" "foodt-1 ties INT" "xit-1 times int"))




*************************************************************************************************************************************************************
*************************************************************************************************************************************************************
*********************************************************Figure 5********************************************************************************************
*************************************************************************************************************************************************************
*************************************************************************************************************************************************************
xtset pairid year
fvset base 2010 year 
ppmlhdfe application_IDP c.foodcrisis#i.year#c.INT 2011bn.year#c.INT 2012bn.year#c.INT 2013bn.year#c.INT 2014bn.year#c.INT 2015bn.year#c.INT 2016bn.year#c.INT 2017bn.year#c.INT 2018bn.year#c.INT $control4nc, absorb( scrid#year destid#year pairid )  d
margins, dydx(INT ) at(year=(2010(1)2018) foodcrisis=(0 1)) predict(xb) noestimcheck
marginsplot, xlabel(,ang(45)) legend(size(*0.90) nobox order(1 "No food crisis" 2 "With food crisis" ) position(6)) graphregion(color(white) lwidth(large)) title("") xtitle("")
 	graph save "Graph" "Figure5.gph", replace
graph export "Figure5.jpg", as(jpg) name("Graph") quality(100) replace


*************************************************************************************************************************************************************
*************************************************************************************************************************************************************
*********************************************************Table 2*********************************************************************************************
*************************************************************************************************************************************************************
*************************************************************************************************************************************************************

local mylistd   application_IDP 
foreach d of local mylistd {

eststo, ti("Level of food insecurity"): ppmlhdfe `d' intinsecure1_orig intinsecure2_orig intinsecure3_orig $control4nc    INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store `d'in1

eststo, ti("Persistence of food insecurity"): ppmlhdfe `d' intFoodInsecure1_orig intFoodInsecure2_orig intFoodInsecure3_orig $control4nc  INTL_BRDR_2011- INTL_BRDR_2018 , absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store `d'in2


}


esttab application_IDPin1	application_IDPin2 using "Table2.rtf", pr2 ///
replace nocons star(* 0.10 ** 0.05 *** 0.01) b(3) se(3)  ///
keep(intinsecure1_orig intinsecure2_orig intinsecure3_orig intFoodInsecure1_orig intFoodInsecure2_orig intFoodInsecure3_orig $control4nc)  ///
order(intinsecure1_orig intinsecure2_orig intinsecure3_orig intFoodInsecure1_orig intFoodInsecure2_orig intFoodInsecure3_orig $control4nc)   ///
label nonotes  s(num_n FE3 FE5 FE6 FE7 IDP fo xt, label("Observations" "Country pair FE" "Origin-year FE" "Destination-year FE" "Border-year FE" "IDP" "foodt-1 ties INT" "xit-1 times int"))

esttab application_IDPin1	application_IDPin2 using "Table2.tex", pr2 ///
replace nocons star(* 0.10 ** 0.05 *** 0.01) b(3) se(3)  ///
keep(intinsecure1_orig intinsecure2_orig intinsecure3_orig intFoodInsecure1_orig intFoodInsecure2_orig intFoodInsecure3_orig $control4nc)  ///
order(intinsecure1_orig intinsecure2_orig intinsecure3_orig intFoodInsecure1_orig intFoodInsecure2_orig intFoodInsecure3_orig $control4nc)   ///
label nonotes  s(num_n FE3 FE5 FE6 FE7 IDP fo xt, label("Observations" "Country pair FE" "Origin-year FE" "Destination-year FE" "Border-year FE" "IDP" "foodt-1 ties INT" "xit-1 times int"))





*************************************************************************************************************************************************************
*************************************************************************************************************************************************************
*********************************************************Figure 7********************************************************************************************
*************************************************************************************************************************************************************
*************************************************************************************************************************************************************


coefplot (application_IDP4,  recast(bar) color(blue) fintensity(40) offset(0))  (application_IDPin1,  recast(bar) color(blue) fintensity(40) offset(0)), ///
      barwidth(0.5)  ///
    keep(intfoodcrisis intinsecure1_orig intinsecure2_orig intinsecure3_orig) vertical ///
    transform(intfoodcrisis intinsecure1_orig intinsecure2_orig intinsecure3_orig = (exp(@)-1)*100) ///
    noci ///
    yline(0) ytitle("Percentage change", size(vsmall)) ///
    ylabel(, labsize(vsmall)) xlabel(, labsize(vsmall)) ///
    ylabel(-50(50)650, labsize(vsmall) nogrid) ///
    graphregion(color(white)) ///
	xline(1.5 10, lp(shortdash) lc(gs6)) ///
    plotregion(margin(b = 0)) ///
    aspect(0.45) nokey  
	*addplot: , legend(order(1 "Ultimate" 2 "Direct") on)  norescaling
 	graph save "Graph" "Figure6.gph", replace
	graph export "Figure6.jpg", as(jpg) name("Graph") quality(100) replace




*************************************************************************************************************************************************************
*************************************************************************************************************************************************************
*********************************************************Table 3*********************************************************************************************
******************************************************Tables A8, A9 & A10************************************************************************************
*************************************************************************************************************************************************************

** Table 3 summarizes the following estimates that aproximate the effect food crises when the destination country is developed or neighbour, or when the source country has a negative GDP pc growth. Full estimates are reported on Tables A8, A9 and A10. The following code resplicates Tables A8, A9 and A10.

/*
local mylistd D_intinsecure1_orig D_intinsecure2_orig D_intinsecure3_orig D_intFoodInsecure1_orig D_intFoodInsecure2_orig D_intpolinstab_only_o D_intmigration_only_o D_intweatherdis_only_o D_intmoreonecause D_intFoodInsecure3_orig D_intfoodcrisis D_inttwocauses D_intmoretntwo D_inteconomic_only_o D_inteco_poli D_inteco_mig D_inteco_weather D_intpoli_mig D_intpoli_weather D_intmig_weather
foreach d of local mylistd {
label var `d' "x Developed"
}
*/

***************************
************ Table A8******
***************************
******Developed************
***************************

local mylistd   application_IDP 
foreach d of local mylistd {

eststo, ti("`d'"): ppmlhdfe `d' intfoodcrisis D_intfoodcrisis $control4nc  INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store `d'd1


eststo, ti("`d'"): ppmlhdfe `d' intinsecure1_orig D_intinsecure1_orig intinsecure2_orig D_intinsecure2_orig intinsecure3_orig D_intinsecure3_orig $control4nc INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store `d'd2

eststo, ti("`d'"): ppmlhdfe `d'  intFoodInsecure1_orig D_intFoodInsecure1_orig intFoodInsecure2_orig D_intFoodInsecure2_orig intFoodInsecure3_orig D_intFoodInsecure3_orig $control4nc INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store `d'd3


}



esttab application_IDPd1	application_IDPd2	application_IDPd3 using "TableA8.rtf", pr2 ///
replace nocons star(* 0.10 ** 0.05 *** 0.01) b(3) se(3)  ///
keep(intfoodcrisis D_intfoodcrisis intinsecure1_orig D_intinsecure1_orig intinsecure2_orig D_intinsecure2_orig intinsecure3_orig D_intinsecure3_orig intFoodInsecure1_orig D_intFoodInsecure1_orig intFoodInsecure2_orig D_intFoodInsecure2_orig intFoodInsecure3_orig D_intFoodInsecure3_orig $control4nc)  ///
order(intfoodcrisis D_intfoodcrisis intinsecure1_orig D_intinsecure1_orig intinsecure2_orig D_intinsecure2_orig intinsecure3_orig D_intinsecure3_orig intFoodInsecure1_orig D_intFoodInsecure1_orig intFoodInsecure2_orig D_intFoodInsecure2_orig intFoodInsecure3_orig D_intFoodInsecure3_orig $control4nc)   ///
label nonotes s(num_n FE3 FE5 FE6 FE7 IDP fo xt, label("Observations" "Country pair FE" "Origin-year FE" "Destination-year FE" "Border-year FE" "IDP" "foodt-1 ties INT" "xit-1 times int"))

esttab application_IDPd1	application_IDPd2	application_IDPd3 using "TableA8.tex", pr2 ///
replace nocons star(* 0.10 ** 0.05 *** 0.01) b(3) se(3)  ///
keep(intfoodcrisis D_intfoodcrisis intinsecure1_orig D_intinsecure1_orig intinsecure2_orig D_intinsecure2_orig intinsecure3_orig D_intinsecure3_orig intFoodInsecure1_orig D_intFoodInsecure1_orig intFoodInsecure2_orig D_intFoodInsecure2_orig intFoodInsecure3_orig D_intFoodInsecure3_orig $control4nc)  ///
order(intfoodcrisis D_intfoodcrisis intinsecure1_orig D_intinsecure1_orig intinsecure2_orig D_intinsecure2_orig intinsecure3_orig D_intinsecure3_orig intFoodInsecure1_orig D_intFoodInsecure1_orig intFoodInsecure2_orig D_intFoodInsecure2_orig intFoodInsecure3_orig D_intFoodInsecure3_orig $control4nc)   ///
label nonotes s(num_n FE3 FE5 FE6 FE7 IDP fo xt, label("Observations" "Country pair FE" "Origin-year FE" "Destination-year FE" "Border-year FE" "IDP" "foodt-1 ties INT" "xit-1 times int"))

   

/*
local mylistd intfoodcrisis_cont intinsecure1_orig_cont intinsecure2_orig_cont intinsecure3_orig_cont intFoodInsecure1_orig_cont intFoodInsecure2_orig_cont intFoodInsecure3_orig_cont
foreach d of local mylistd {
label var `d' "x Contiguity"
}
*/
***************************
************ Table A9******
***************************
****Neighbour**************
***************************

eststo, ti("`d'"): ppmlhdfe application_IDP intfoodcrisis intfoodcrisis_cont $control4nc INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store cont1

eststo, ti("`d'"): ppmlhdfe application_IDP intinsecure1_orig intinsecure1_orig_cont intinsecure2_orig intinsecure2_orig_cont intinsecure3_orig intinsecure3_orig_cont  $control4nc INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store cont2

eststo, ti("`d'"): ppmlhdfe application_IDP intFoodInsecure1_orig intFoodInsecure1_orig_cont intFoodInsecure2_orig intFoodInsecure2_orig_cont  intFoodInsecure3_orig intFoodInsecure3_orig_cont  $control4nc INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store cont3



esttab cont1 cont2 cont3 using "TableA9.rtf", pr2 ///
replace nocons star(* 0.10 ** 0.05 *** 0.01) b(3) se(3)  ///
keep(intfoodcrisis intfoodcrisis_cont intinsecure1_orig intinsecure1_orig_cont intinsecure2_orig intinsecure2_orig_cont intinsecure3_orig intinsecure3_orig_cont  intFoodInsecure1_orig intFoodInsecure1_orig_cont intFoodInsecure2_orig intFoodInsecure2_orig_cont  intFoodInsecure3_orig intFoodInsecure3_orig_cont $control4nc)  ///
order(intfoodcrisis intfoodcrisis_cont intinsecure1_orig intinsecure1_orig_cont intinsecure2_orig intinsecure2_orig_cont intinsecure3_orig intinsecure3_orig_cont  intFoodInsecure1_orig intFoodInsecure1_orig_cont intFoodInsecure2_orig intFoodInsecure2_orig_cont  intFoodInsecure3_orig intFoodInsecure3_orig_cont $control4nc)   ///
label nonotes s(num_n FE3 FE5 FE6 FE7 IDP fo xt, label("Observations" "Country pair FE" "Origin-year FE" "Destination-year FE" "Border-year FE" "IDP" "foodt-1 ties INT" "xit-1 times int"))

esttab cont1 cont2 cont3 using "TableA9.tex", pr2 ///
replace nocons star(* 0.10 ** 0.05 *** 0.01) b(3) se(3)  ///
keep(intfoodcrisis intfoodcrisis_cont intinsecure1_orig intinsecure1_orig_cont intinsecure2_orig intinsecure2_orig_cont intinsecure3_orig intinsecure3_orig_cont  intFoodInsecure1_orig intFoodInsecure1_orig_cont intFoodInsecure2_orig intFoodInsecure2_orig_cont  intFoodInsecure3_orig intFoodInsecure3_orig_cont $control4nc)  ///
order(intfoodcrisis intfoodcrisis_cont intinsecure1_orig intinsecure1_orig_cont intinsecure2_orig intinsecure2_orig_cont intinsecure3_orig intinsecure3_orig_cont  intFoodInsecure1_orig intFoodInsecure1_orig_cont intFoodInsecure2_orig intFoodInsecure2_orig_cont  intFoodInsecure3_orig intFoodInsecure3_orig_cont $control4nc)   ///
label nonotes s(num_n FE3 FE5 FE6 FE7 IDP fo xt, label("Observations" "Country pair FE" "Origin-year FE" "Destination-year FE" "Border-year FE" "IDP" "foodt-1 ties INT" "xit-1 times int"))


***************************
************ Table A10*****
***************************
****** Growth**************
***************************

eststo, ti("`d'"): ppmlhdfe application_IDP intfoodcrisis intfoodcrisis_growth $control4nc INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store grw1

eststo, ti("`d'"): ppmlhdfe application_IDP intinsecure1_orig intinsecure1_orig_growth intinsecure2_orig intinsecure2_orig_growth intinsecure3_orig intinsecure3_orig_growth  $control4nc INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store grw2

eststo, ti("`d'"): ppmlhdfe application_IDP  intFoodInsecure1_orig intFoodInsecure1_orig_growth intFoodInsecure2_orig intFoodInsecure2_orig_growth  intFoodInsecure3_orig intFoodInsecure3_orig_growth    $control4nc INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store grw3



esttab grw1 grw2 grw3 using "TableA10.rtf", pr2 ///
replace nocons star(* 0.10 ** 0.05 *** 0.01) b(3) se(3)  ///
keep(intfoodcrisis intfoodcrisis_growth intinsecure1_orig intinsecure1_orig_growth intinsecure2_orig intinsecure2_orig_growth intinsecure3_orig intinsecure3_orig_growth  intFoodInsecure1_orig intFoodInsecure1_orig_growth intFoodInsecure2_orig intFoodInsecure2_orig_growth  intFoodInsecure3_orig intFoodInsecure3_orig_growth $control4nc)  ///
order( intfoodcrisis intfoodcrisis_growth intinsecure1_orig intinsecure1_orig_growth intinsecure2_orig intinsecure2_orig_growth intinsecure3_orig intinsecure3_orig_growth  intFoodInsecure1_orig intFoodInsecure1_orig_growth intFoodInsecure2_orig intFoodInsecure2_orig_growth  intFoodInsecure3_orig intFoodInsecure3_orig_growth $control4nc)   ///
label nonotes s(num_n FE3 FE5 FE6 FE7 IDP fo xt, label("Observations" "Country pair FE" "Origin-year FE" "Destination-year FE" "Border-year FE" "IDP" "foodt-1 ties INT" "xit-1 times int"))

esttab grw1 grw2 grw3 using "TableA10.tex", pr2 ///
replace nocons star(* 0.10 ** 0.05 *** 0.01) b(3) se(3)  ///
keep(intfoodcrisis intfoodcrisis_growth intinsecure1_orig intinsecure1_orig_growth intinsecure2_orig intinsecure2_orig_growth intinsecure3_orig intinsecure3_orig_growth  intFoodInsecure1_orig intFoodInsecure1_orig_growth intFoodInsecure2_orig intFoodInsecure2_orig_growth  intFoodInsecure3_orig intFoodInsecure3_orig_growth $control4nc)  ///
order( intfoodcrisis intfoodcrisis_growth intinsecure1_orig intinsecure1_orig_growth intinsecure2_orig intinsecure2_orig_growth intinsecure3_orig intinsecure3_orig_growth  intFoodInsecure1_orig intFoodInsecure1_orig_growth intFoodInsecure2_orig intFoodInsecure2_orig_growth  intFoodInsecure3_orig intFoodInsecure3_orig_growth $control4nc)   ///
label nonotes s(num_n FE3 FE5 FE6 FE7 IDP fo xt, label("Observations" "Country pair FE" "Origin-year FE" "Destination-year FE" "Border-year FE" "IDP" "foodt-1 ties INT" "xit-1 times int"))



****************************************************************************************************************************************************************
****************************************************************************************************************************************************************
***** Tables in the appendix************************************************************************************************************************************
****************************************************************************************************************************************************************
****************************************************************************************************************************************************************

*************************************************************************************************************************************************************
*************************************************************************************************************************************************************
*********************************************************Table A4********************************************************************************************
*************************************************************************************************************************************************************
*************************************************************************************************************************************************************

eststo, ti("`d'"): ppmlhdfe application_IDP  intfoodcrisis $control4nc2 $bilatt INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local ofe ""
estadd local dfe ""
estadd local pfe "X"
estadd local yrfe ""
estadd local oyfe "X"
estadd local dyfe "X"
estadd local brd "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store s1

eststo, ti("`d'"): ppmlhdfe application_IDP intfoodcrisis $control4nc2 $bilatt, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local ofe ""
estadd local dfe ""
estadd local pfe "X"
estadd local yrfe ""
estadd local oyfe "X"
estadd local dyfe "X"
estadd local brd ""
estadd local fo "Yes"
estadd local xt "Yes"
est store s2

eststo, ti("`d'"): ppmlhdfe application_IDP intfoodcrisis $control4nc2 $bilatt $bilat, absorb( scrid#year destid#year) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local ofe ""
estadd local dfe ""
estadd local pfe ""
estadd local yrfe ""
estadd local oyfe "X"
estadd local dyfe "X"
estadd local brd ""
estadd local fo "Yes"
estadd local xt "Yes"
est store s3

eststo, ti("`d'"): ppmlhdfe application_IDP  $bilatt $bilat, absorb( scrid#year destid#year) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local ofe ""
estadd local dfe ""
estadd local pfe ""
estadd local yrfe ""
estadd local oyfe "X"
estadd local dyfe "X"
estadd local brd ""
estadd local fo "No"
estadd local xt "No"
est store s4


eststo, ti("`d'"): ppmlhdfe application_IDP foodcrisis $country $control4ncnoint $bilatt $bilat, absorb( scrid destid year) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local ofe "X"
estadd local dfe "X"
estadd local pfe ""
estadd local yrfe "X"
estadd local oyfe ""
estadd local dyfe ""
estadd local brd ""
estadd local fo "No"
estadd local xt "No"
est store s5





eststo, ti("`d'"): ppmlhdfe application_IDP foodcrisis $country $control4ncnoint $bilatt $bilat if  ccode_dest!= ccode_scr, absorb(destid#year pairid) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "No"
estadd local ofe ""
estadd local dfe ""
estadd local pfe "X"
estadd local yrfe ""
estadd local oyfe ""
estadd local dyfe "X"
estadd local brd ""
estadd local fo "No"
estadd local xt "No"
est store s6

eststo, ti("`d'"): ppmlhdfe application_IDP foodcrisis $country $control4ncnoint $bilatt $bilat if  ccode_dest!= ccode_scr, absorb(pairid year) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "No"
estadd local ofe ""
estadd local dfe ""
estadd local pfe "X"
estadd local yrfe "X"
estadd local oyfe ""
estadd local dyfe ""
estadd local brd ""
estadd local fo "No"
estadd local xt "No"
est store s7


eststo, ti("`d'"): ppmlhdfe application_IDP foodcrisis $country $control4ncnoint $bilatt $bilat if  ccode_dest!= ccode_scr, absorb(scrid destid year) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "No"
estadd local ofe "X"
estadd local dfe "X"
estadd local pfe ""
estadd local yrfe "X"
estadd local oyfe ""
estadd local dyfe ""
estadd local brd ""
estadd local fo "No"
estadd local xt "No"
est store s8


esttab s1 s2 s3 s4 s5 s6 s7 s8 using "TableA4.rtf", pr2 ///
replace nocons star(* 0.10 ** 0.05 *** 0.01) b(3) se(3)  ///
keep(intfoodcrisis foodcrisis $control4nc2 $control4ncnoint $country $bilatt $bilat)  ///
order(intfoodcrisis foodcrisis $control4nc2 $control4ncnoint $country $bilatt $bilat)   ///
label nonotes s(num_n ofe dfe pfe yrfe oyfe dyfe brd IDP fo xt, label("Observations" "Origin FE" "Destination FE" "Country pair FE" "Year FE"  "Origin-year FE" "Destination-year FE" "Border-year FE" "IDP" "foodt-1 ties INT" "xit-1 times int"))


esttab s1 s2 s3 s4 s5 s6 s7 s8 using "TableA4.tex", pr2 ///
replace nocons star(* 0.10 ** 0.05 *** 0.01) b(3) se(3)  ///
keep(intfoodcrisis foodcrisis $control4nc2 $control4ncnoint $country $bilatt $bilat)  ///
order(intfoodcrisis foodcrisis $control4nc2 $control4ncnoint $country $bilatt $bilat)   ///
label nonotes s(num_n ofe dfe pfe yrfe oyfe dyfe brd IDP fo xt, label("Observations" "Origin FE" "Destination FE" "Country pair FE" "Year FE"  "Origin-year FE" "Destination-year FE" "Border-year FE" "IDP" "foodt-1 ties INT" "xit-1 times int"))






*************************************************************************************************************************************************************
*************************************************************************************************************************************************************
*********************************************************Table A5********************************************************************************************
*************************************************************************************************************************************************************
*************************************************************************************************************************************************************


eststo, ti("application_IDP"): ppmlhdfe application_IDP  $control4 intdisaster intL1polstab intL1unemploy_s INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store application_IDP5

esttab application_IDP5	 using "TableA5.tex", pr2 ///
replace nocons star(* 0.10 ** 0.05 *** 0.01) b(3) se(3)  ///
keep($controltab intdisaster intL1polstab intL1unemploy_s)  ///
order($controltab intdisaster intL1polstab intL1unemploy_s)   ///
label nonotes s(num_n FE3 FE5 FE6 FE7 IDP fo xt, label("Observations" "Country pair FE" "Origin-year FE" "Destination-year FE" "Border-year FE" "IDP" "foodt-1 ties INT" "xit-1 times int"))


esttab application_IDP5	 using "TableA5.rtf", pr2 ///
replace nocons star(* 0.10 ** 0.05 *** 0.01) b(3) se(3)  ///
keep($controltab intdisaster intL1polstab intL1unemploy_s)  ///
order($controltab intdisaster intL1polstab intL1unemploy_s)   ///
label nonotes s(num_n FE3 FE5 FE6 FE7 IDP fo xt, label("Observations" "Country pair FE" "Origin-year FE" "Destination-year FE" "Border-year FE" "IDP" "foodt-1 ties INT" "xit-1 times int"))


*************************************************************************************************************************************************************
*************************************************************************************************************************************************************
*********************************************************Table A6********************************************************************************************
*************************************************************************************************************************************************************
*************************************************************************************************************************************************************

eststo, ti("`d'"): ppmlhdfe application_IDP  $control1q  INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store qf1



eststo, ti("`d'"): ppmlhdfe application_IDP  $control2q INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store qf2



eststo, ti("`d'"): ppmlhdfe application_IDP  $control3q  INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store qf3


eststo, ti("`d'"): ppmlhdfe application_IDP  $control4q  INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store qf4



esttab qf1 qf2 qf3 qf4	 using "TableA6.tex", pr2 ///
replace nocons star(* 0.10 ** 0.05 *** 0.01) b(3) se(3)  ///
keep($controltabq)  ///
order($controltabq)   ///
label nonotes s(num_n FE3 FE5 FE6 FE7 IDP fo xt, label("Observations" "Country pair FE" "Origin-year FE" "Destination-year FE" "Border-year FE" "IDP" "foodt-1 ties INT" "xit-1 times int"))

esttab qf1 qf2 qf3 qf4	 using "TableA6.rtf", pr2 ///
replace nocons star(* 0.10 ** 0.05 *** 0.01) b(3) se(3)  ///
keep($controltabq)  ///
order($controltabq)   ///
label nonotes s(num_n FE3 FE5 FE6 FE7 IDP fo xt, label("Observations" "Country pair FE" "Origin-year FE" "Destination-year FE" "Border-year FE" "IDP" "foodt-1 ties INT" "xit-1 times int"))



*************************************************************************************************************************************************************
*************************************************************************************************************************************************************
*********************************************************Table A7********************************************************************************************
*************************************************************************************************************************************************************
*************************************************************************************************************************************************************



eststo, ti("`d'"): ppmlhdfe application_IDP  L2intfoodcrisis $controltab2  INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store lag4

eststo, ti("`d'"): ppmlhdfe application_IDP  L2intfoodcrisis INTL_BRDR_2011- INTL_BRDR_2018 if _est_lag4==1, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store lag1



eststo, ti("`d'"): ppmlhdfe application_IDP  L2intfoodcrisis $controllag1 INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store lag2



eststo, ti("`d'"): ppmlhdfe application_IDP  L2intfoodcrisis  $controlag2  INTL_BRDR_2011- INTL_BRDR_2018, absorb( scrid#year destid#year pairid ) vce(cluster pairid)
local k=e(N)
local r=  e(r2_p)
estadd local num_n `k'
estadd local r2 `r'
estadd local IDP "Yes"
estadd local FE1 ""
estadd local FE2 ""
estadd local FE3 "X"
estadd local FE4 ""
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
estadd local fo "Yes"
estadd local xt "Yes"
est store lag3





esttab lag1 lag2 lag3 lag4	 using "TableA7.tex", pr2 ///
replace nocons star(* 0.10 ** 0.05 *** 0.01) b(3) se(3)  ///
keep(L2intfoodcrisis $controltab2)  ///
order(L2intfoodcrisis $controltab2)   ///
label nonotes s(num_n FE3 FE5 FE6 FE7 IDP fo xt, label("Observations" "Country pair FE" "Origin-year FE" "Destination-year FE" "Border-year FE" "IDP" "foodt-2 ties INT" "xit-2 times int"))

esttab lag1 lag2 lag3 lag4	 using "TableA7.rtf", pr2 ///
replace nocons star(* 0.10 ** 0.05 *** 0.01) b(3) se(3)  ///
keep(L2intfoodcrisis $controltab2)  ///
order(L2intfoodcrisis $controltab2)   ///
label nonotes s(num_n FE3 FE5 FE6 FE7 IDP fo xt, label("Observations" "Country pair FE" "Origin-year FE" "Destination-year FE" "Border-year FE" "IDP" "foodt-2 ties INT" "xit-2 times int"))


*************************************************************************************************************************************************************
*************************************************************************************************************************************************************
*********************************************************Table 4*********************************************************************************************
**********************************************************IV*************************************************************************************************
*************************************************************************************************************************************************************


** Instrument
gen i_instrument_USwheatstock=INT*instrument_USwheatstock


*****First step regression (for relevance of instrument)
regress foodcrisis instrument_USwheatstock i.ISO3_orig i.year


********   Predicting fixed effects derived from PPML for second stage *******

ppml_panel_sg application_IDP rta lnmigration if  ccode_scr!=ccode_dest, ex(scrid) imp(destid) year(year) genM(e_dest_ext_FE_d) genO(scr_yr_int)
local k=e(N)
capture local r=  e(r2_p)
estadd local num_n `k'
capture estadd local r2 `r'
estadd local FE5 "X"
estadd local FE6 "X"
estadd local FE7 "X"
est store base


duplicates drop scr_yr_int, force

******* Estimations with US_wheat stock and instrument ***

global controls4 intany_crisis intL1VoiceandAccountability intL1death_pop intgdppcgr 


reghdfe scr_yr_int intfoodcrisis, absorb (scrid year) cluster(scryr)
local k=e(N)
estadd local num_n `k'
estadd local FE1 "X"
estadd local FE2 "X"
estadd local dep "$\\hat{\\lambda}_{it}$"
*estadd local T1 ""
*estadd local T2 ""
est store re2


ivreghdfe scr_yr_int  (intfoodcrisis=instrument_USwheatstock), first savefprefix(firstiv) absorb(scrid year)  cluster(scryr)
local k=e(N)
estadd local num_n `k'
estadd local FE1 "X"
estadd local FE2 "X"
estadd scalar T1 `e(j)'
estadd scalar T2 `e(j)'
estadd scalar T3 `e(idstat)'
estadd scalar T4 `e(idp)'
*estadd local T1 `e(cdf)'
*estadd local T2 `e(sstat)'
estadd scalar T5 `e(F)'
estadd scalar T6 `e(Fp)'
estadd local dep "$\\hat{\\lambda}_{it}$"
est store ivf1
estimates restore firstivintfoodcrisis
estimates store firstivintfoodcrisis
**OJO, al parecer el comando no es capaz de guardar este local tras recuperar los resultados del primer stage.
*Hay que poner la dependent variable a mano en la tabla. 
estadd local dep "Food crisis"



ivreghdfe scr_yr_int $controls4  (intfoodcrisis=instrument_USwheatstock), savefprefix(firstiv2) absorb(scrid year) cluster(scryr)
local k=e(N)
estadd local num_n `k'
estadd local FE1 "X"
estadd local FE2 "X"
estadd scalar T1 `e(j)'
estadd scalar T2 `e(j)'
estadd local dep "$\\hat{\\lambda}_{it}$"
*estadd local T1 `e(cdf)'
*estadd local T2 `e(sstat)'
**OJO, los dos siguientes son stats de la first stage, pero una vez restauradas las estimaciones de la primera
*etapa no permite almacenarlas. Cuidado porque en la tabla se despliegan en la segunda etapa y hay
*que moverlas a mano.
estadd scalar T3 `e(idstat)'
estadd scalar T4 `e(idp)'
est store ivf3
estimates restore firstiv2intfoodcrisis
estimates store firstiv2intfoodcrisis
estadd local dep "$\gamma FOOD_{it-1}$"




reghdfe scr_yr_int  intfoodcrisis $controls4, absorb (scrid year) cluster(scryr)
local k=e(N)
estadd local num_n `k'
estadd local FE1 "X"
estadd local FE2 "X"
estadd local T1 ""
estadd local T2 ""
estadd local T3 ""
estadd local T4 ""
estadd local dep "$\\hat{\\lambda}_{it}$"
est store re9

reghdfe scr_yr_int instrument_USwheatstock, absorb (scrid year) cluster(scryr)
local k=e(N)
estadd local num_n `k'
estadd local FE1 "X"
estadd local FE2 "X"
*estadd local T1 ""
*estadd local T2 ""
estadd local dep "$\\hat{\\lambda}_{it}$"
est store reduced



cd "C:\Users\Usuario\OneDrive\AA Research\Working research projects\Migration and food crisis\Submission\Replication"
esttab reduced firstivintfoodcrisis ivf1 re2 ivf3 re9  using "Table4.rtf", pr2 ///
replace nocons mti("Reduced form" "First stage IV" "Second stage IV" "Second stage no IV" "Second stage IV + controls" "Second stage no IV + controls") star(* 0.10 ** 0.05 *** 0.01) b(3) se(3) ///
keep(instrument_USwheatstock intfoodcrisis $controls4)  ///
order(instrument_USwheatstock intfoodcrisis $controls4)   ///
label nonotes s(num_n FE1 FE2 FE5 FE6 FE7 dep T3 T4, fmt(%9.3f) label("Observations" "Country FEs" "Year FEs" "Origin-year FEs" "Destination-year FEs" "Pair FEs" "Dependent variable" "Kleibergen-Paap F-test" "p-value" ))

esttab reduced firstivintfoodcrisis ivf1 re2 ivf3 re9  using "Table4.tex", pr2 ///
replace nocons mti("Reduced form" "First stage IV" "Second stage IV" "Second stage no IV" "Second stage IV + controls" "Second stage no IV + controls") star(* 0.10 ** 0.05 *** 0.01) b(3) se(3) ///
keep(instrument_USwheatstock intfoodcrisis $controls4)  ///
order(instrument_USwheatstock intfoodcrisis $controls4 )   ///
label nonotes s(num_n FE1 FE2 FE5 FE6 FE7 dep T3 T4 , fmt(%9.3f)  label("Observations" "Country FEs" "Year FEs" "Origin-year FEs" "Destination-year FEs" "Pair FEs" "Dependent variable" "Kleibergen-Paap F-test" "p-value" ))

